Universe With Time Dependent Deceleration Parameter 
and A Term in General Relativity 



Anirudh Pradhan 1 

Department of Mathematics, Hindu Post-graduate College, Zamania-232 331, 

Ghazipur, India 

' E-mail: pradhan@iucaa.emet. in, acpradhan@yahoo.com 

o ' 
o 



00 
(N 



o 
o 



- 1-^ 

x 



Saeed Otarod 

Department of Physics, Yasouj University, Yasouj, Iran 
E-mail: sotarod@mail.yu.ac.ir, sotarod@yahoo.com 



■ Abstract 

CI A new class of exact solutions of Einstein's field equations with perfect 

fluid for an LRS Bianchi type-I spacetime is obtained by using a time 
dependent deceleration parameter. We have obtained a general solution 
of the field equations from which three models of the universe are derived: 
exponential, polynomial and sinusoidal form respectively. The behaviour 
of these models of the universe are also discussed in the frame of reference 
^ of recent supernovae la observations. 
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1 Introduction 

The Bianchi cosmologies play an important role in theoretical cosmology and 
have been much studied since the 1960s. A Bianchi cosmology represents a spa- 
tially homogeneous universe, since by definition the spacetime admits a three- 
parameter group of isometries whose orbits are spacelike hyper-surfaces. These 
models can be used to analyze aspects of the physical Universe which pertain 
to or which may be affected by anisotropy in the rate of expansion, for exam- 
ple , the cosmic microwave background radiation, nucleosynthesis in the early 
universe, and the question of the isotropization of the universe itself £Q. For 
simplification and description of the large scale behaviour of the actual universe, 
locally rotationally symmetric [henceforth referred as LRS] Bianchi I spacetime 
have widely studied — EH- When the Bianchi I spacetime expands equally in 
two spatial directions it is called locally rotationally symmetric. These kinds of 
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models are interesting because Lidsey [7] showed that they are equivalent to a 
flat (FRW) universe with a self-interesting scalar field and a free massless scalar 
field, but produced no explicit example. Some explicit solutions were pointed 
out in references [HI ■ 

The Einstein's field equations are coupled system of high non-linear differ- 
ential equations and we seek physical solutions to the field equations for their 
applications in cosmology and astrophysics. In order to solve the field equations 
we normally assume a form for the matter content or that spacetime admits 
killing vector symmetries ^UJ- Solutions to the field equations may also be 
generated by applying a law of variation for Hubble's parameter which was pro- 
posed by Berman ^U- hi simple cases the Hubble law yields a constant value 
of deceleration parameter. It is worth observing that most of the well-known 
models of Einstein's theory and Brans-Deke theory with curvature parameter 
k = 0, including inflationary models, are models with constant deceleration pa- 
rameter. In earlier literature cosmological models with a constant deceleration 
parameter have been studied by Berman Berman and Gomide Johri 
and Dcsikan Singh and Desikan Maharaj and Naidoo [I^], Pradhan 
and et al. |16j and others. But redshift magnitude test has had a chequered 
history. During the 1960s and the 1970s, it was used to draw very categorical 
conclusions. The deceleration parameter q was then claimed to lie between 
and 1 and thus it was claimed that the universe is decelerating. Today's situ- 
ation, we feel, is hardly different. Observations ^1E3 of Type la Supernovae 
(SNe) allow to probe the expansion history of the universe. The main conclusion 
of these observations is that the expansion of the universe is accelerating. So 
we can consider the cosmological models with variable deceleration parameter. 
The readers are advised to see the papers by Vishwakarma and Narlikar [T§| 
and Virey et al. |20j and references therein for a review on the determination 
of the deceleration parameter from Supernovae data. 

Motivated with the situation discussed above, in this paper we can focus 
upon the problem of establishing a formalism for studying the relativistic evo- 
lution for a time dependent deceleration parameter in an expanding universe. 
This paper is organized as follows. The metric and the field equations are pre- 
sented in Section 2. In Section 3 we deal with a general solution. The Sections 
4, 5, and 6 contain the three different cases for the solutions in exponential, 
polynomial and sinusoidal forms respectively. Finally in Section 7 concluding 
remarks will be given. 

2 The Metric and Field Equations 

We consider the LRS Bianchi type-I metric in the form [2] 

ds 2 = dt 2 -A 2 dx 2 -B 2 {dy 2 + dz 2 ), (1) 
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where A and B are functions of x and t. The energy momentum-tensor in the 
presence of perfect fluid has the form 

Tij = (p + p)uiUj - pg lj , (2) 

where p, p are the energy density thermodynamical pressure respectively and 
m is the four velocity vector satisfying the relations 

u t u l = 1, (3) 

The Einstein's field equations (in gravitational units c = 1, G = 1) read as 

R io ~ \ R 9ij + Aftj = SnTij, (4) 

where Rij is the Ricci tensor; R = g y Rij is the Ricci scalar. The Einstein's 
field equations Q for the line element (|TJ has been set up as 
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The energy conservation equation yields 



P+ip + p) (A + lf) + A = °. ( 9 ) 

where dots and primes indicate partial differentiation with respect to t and x 
respectively. 

In order to completely determine the system, we choose a barotropic equation 
of state 

P = 7P,0<7<1. (10) 

3 Solution of the Field Equations 

Equation |JBJ, after integration, yields 

A = ^r, (ii) 
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where I is an arbitrary function of x. Equations (JSJ and (J7J, with the use of Eq. 
l(TT|) . reduces to 



B d B 



B d (B 



B' dx \ B 
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(12) 



If we assume ^- to be a function of x alone, then A and B are separable in x 
and t. Hence, after integrating Eq. 1|12|) gives 



B = tS{t), 



(13) 



where S is a scale factor which is an arbitrary function of t. Thus from Eqs. 
l(TT|) and ljT3jl . we have 

4 = j& (14) 
Now the metric Q is reduced to the form 

ds 2 = dt 2 - S 2 [dX 2 + e 2X (dy 2 + dz 2 )} , (15) 
where X = In £. The mass-density, pressure and Ricci scalar are obtained as 



8np ^s 2 
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(16) 

(17) 
(18) 



The function S(t) remains undetermined. To obtain its explicit dependence on 
t, one may have to introduce additional assumption. To achieve this, we assume 
the deceleration parameter to be variable, i.e. 



1p 



H + H 2 
H 2 



b (variable), 



(19) 



where H = ^ is the Hubble parameter. The above equation may be rewritten 



as 



S ,S 2 



s s 2 

The general solution of Eq. 12U|) is given by 

J ^ dS dS = 



0. 



dS = t + k, 



(20) 



(21) 



where k is an integrating constant. 

In order to solve the problem completely, we have to choose J ^dS in such 
a manner so that Eq. (|21|l be integrable. 
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Let us consider 

J^dS = lnL(S), (22) 

which does not effect the nature of generality of solution. Hence from Eqs. Ij21|l 
and (J22J, 

one can obtain 

J L(S)dS = t + k. (23) 

Of course the choice of L(S) is quite arbitrary but, since we are looking for phys- 
ically viable models of the universe consistent with observations, we consider the 
following three cases: 



4 Solution in the Exponential Form 

Let us consider L(S) = -j^g, where k\ is an arbitrary constant. 
In this case, on integrating, Eq. (|23[) gives the exact solution 

S(t) = k 2 e klt , (24) 

where k 2 is an arbitrary constant. Using Eas. (|ll)fl and 124fl in Eqs. Hltj[) - (|18[) . 

the mass-density, cosmological term and Rlcci scalar are obtained as 

rrt _/(37+l) __3_\ _L 2fcf 1 / (3^ + 1) 3) a 
P{ > \ 7 + 1 87rJ fc 2 2 e 2fc i* (l+ 7 )fc 2 e fe i* \ j+1 8tt j u { ' 



A:|e 2fc i* J \l+7/ (l + 7)fc 2 e fel * 

From Eq. I|24|) . since scale factor can not be negative, we find S(t) is positive 
if k% > 0. From Figure ^ it can be deduced that at the early stages of the 
universe i.e. near t = 0, the scale factor of the universe had been approximately 
constant and had increased very slowly. At an specific time the universe has 
exploded suddenly and it has expanded to large scale. This fits nicely with Big 
Bang scenario. 

From Eqs. and l(2*B|l. it is observed that p(t) > and A(t) > for 

0<i<ooif0<fc2<l. Figure [5] clearly shows this behaviour of p(t). 

From Eq. (|26J) . we observe that the cosmological term is a decreasing func- 
tion of time and it approaches a small value as time progresses (i.e. the present 
epoch), which explains the small and positive value of A at present (Pcrlmuttcr 
et al. Riess et al. Garnavich et al. [23]; Schmidt et al. 24 ). Figure 
|3j clearly shows this behaviour of A as decreasing function of time. From Eq. 
Il27|) . we see that the Ricci scalar remain positive for 

ki > 1 =. 

y/k 2 (l + k 2 ) 
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Figure 1: The plot of scale factor S(t) vs time with parameters k\ = 0.01, 
fc 2 = 0.5, and 7 = 0.5 
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Figure 2: The plot of energy density p(t) vs time with parameters k\ = 0.01, 
k 2 = 0.5, and 7 = 0.5 
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Figure 3: The plot of cosmological term (A) vs time with parameters k\ = 0.01, 
k-2 = 0.5, and 7 = 0.5 



5 Solution in the Polynomial Form 

Let L(S) = — — 7=j , where and k^ are constants. 

In this case, on integrating, Eq. (|23[) gives the exact solution 

S(t) = ai t 2 + a 2 t + a 3 , (28) 

where oti, a 2 and a 3 are arbitrary constants. Using Eqs. (|10fl and (|28|> in Eqs. 
I|l(j|) - (|18fl . the mass-density, cosmological term and Ricci scalar are obtained as 



_ [(l + 3 7 )-(l + 37)(2a 1 t + a 2 ) 2 -4a 1J 
P[) ~ (l + 7 )(aii 2 +a 2 t + a 3 ) 2 



[(2 ttl t + a 2 )='-2a 1 -l] 
HJ 47r(l + 7 )(a 1 i 2 + a 2 i + a 3 ) 2 ' ^ 
R = 6[(2a 1 t + a 2 ) 2 + 2a 1 - 1] 
(ait 2 + a 2 t + a 3 ) 2 

From Eq. (|28|l . it is observed that S(t) > for < t < 00 if a±, a 2 and a 3 
are positive constants. From Figure 4, it is observed that the scale factor is a 
decreasing function of time which means that our universe is expanding. 

In order to have p > for < t < 00, we must have a 2 > 2a± + 1. From Eq. 
<|31|) we observe that Ricci scalar remains positive if a 2 > 1-201 . This condition 
also implies that a± < \. Figure 5 clearly shows the decreasing behaviour of 
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Figure 4: The plot of scale factor S(t) vs time with parameters a\ = 1.00, 
a 2 = 4.00, a 3 = 1.00 and 7 = 0.5 
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Figure 5: The plot of energy density p(t) vs time with parameters a.\ = 1.00, 
a 2 = 4.00, a 3 = 1.00 and 7 = 0.5 
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Figurc 6: The plot of cosmological term (A) vs time with parameters ot\ = 1.00, 
a 2 = 4.00, a 3 = 1.00 and 7 = 0.5 



p(t) as time increases and is always positive. The interesting point is that all 
physical parameters in our models arc defined at t = and we do not have any 
singularity. 

It is observed from Eq. that A(t) remains always negative but decreasing 
function of time. From the Figure it can be seen the behaviour of A as a 
decreasing function of time. By decreasing we mean its absolute magnitude 
approaches zero which is acceptable physically. A negative cosmological term 
adds to the attractive gravity of matter; therefore, universe with a negative 
cosmological term is invariably doomed to recollapse. A positive cosmological 
term resists the attractive gravity of matter due to its negative pressure. For 
most universe cosmological term eventually dominates over the attraction of 
matter and drives the universe to expands exponentially. 

6 Solution in the Sinusoidal Form 

If we set L(S) = — ^===g , where (3 is constant. 

In this case, on integrating, Eq. (12,'ill gives the exact solution 



S = Msm((3t) + Ncos(fJt) + ft, 



(32) 



infinity 



Figure 7: The plot of scale factor S(t) vs time with parameters M 
N = 1.00, = 10.00, 0i = 0.2, and 7 = 0.5 



2.00, 



where M, N and 0i are constants. Using Eqs.JTJJ and in E qs. l(TB|) -(tT8l). 

the mass-density, cosmological term and Ricci scalar are obtained as 



4tt(1 +7)p = 



2 (Mcos(0t) -Nsin(0t)) 2 + 2 (M sm(0t) - N cos(0it)) - 1 



(Msin(/3i) + Ncos(0t) + 0i 



(33) 



(1+7)A = 



R 



(1 + 3 7 )/3 2 (Mcos(/3t) -Nsm(0t)) 2 - (1 + 3 7 ) - 2k\ (Msin(0t) + N cos(0t)) 



(Msm(0t) +Ncos(0t) + 0iY 
(M0cos(0t) - N0sin{0t)) 2 - 2 {Msin(0t) + N cos(0t)) - 1 



(Msm(0t) +Ncos(0t) + 0iY 



(34) 
(35) 



Since, in this case, we have many alternatives for choosing values of M, N, 0. 
0i, it is just enough to look for suitable values of these parameters, such that 
the physical initial and boundary conditions are satisfied. We are trying to find 
feasible interpretation and situations relevant to this case. Further study in this 
case is in progress. 

From the Figure0it is observed that at early stages of the universe, the scale 
of the universe increases gently and then decreases sharply, and afterwords it 
will oscillate for ever. 

From Figure |H1 and Figure we conclude that at early stages of the universe 
the matter is created as a result of loss of vacuum energy and at a particular 
epoch it has started to oscillate for ever due to sinusoidal property. It is worth 
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Figure 8: The plot ol energy density p(t) vs time with parameters M — 2.00, 
N = 1.00, /3 = 10.00, 0! = 0.2, and 7 = 0.5 
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Figure 9: The plot of cosmological term (A) vs time with parameters M = 2.00, 
N = 1.00, /3 = 10.00, /3i = 0.2, and 7 = 0.5 
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to mention here that in this case oscillation takes place in positive quadrant 
which has physical meaning. 

7 Conclusions 

In this paper we have described a new class of LRS Bianchi type I cosmological 
models with a perfect fluid as the source of matter by applying a variable de- 
celeration parameter. Generally, the models are expanding, non-shearing and 
isotropic in nature. 

The cosmological term in the model in Section 4 is a decreasing function of 
time and this approaches a small value as time increases (i.e. present epoch). 
The value of the cosmological "term" for this model is found to be small and 
positive which is supported by the results from recent supernovae observations 
obtained by the High-Z Supernova Team and Supernova Cosmological Project 
(Perlmutter et al. Riess et al. (221) Garnavich et al. [531; Schmidt et 

al. |24|V The cosmological term in the model in Section 5 is also a decreasing 
function of time but it is always negative. A negative cosmological term adds to 
the attractive gravity of matter; therefore, universe with a negative cosmological 
term is invariably doomed to recollapse. The cosmological term in Section 6 also 
decreases while time increases to a specific instant. During this period as we 
can understand from Figure [51 we will have enough matter creation to force the 
universe to oscillate for ever due to sinusoidal property of A. This means we 
always have annihilation and creation of matter permanently. At this point one 
more sentence may be added to our discussion and i.e. as the graphs for A and 
p in this case the explosion of the universe at the early stages of its creation has 
been only a consequence of matter creation. 
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